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Using high resolution X-Ray diffraction (XRD) on high purity powders, we resolved the structure 
and ab symmetry of the intriguing compound S^VCU from room temperature down to 20 K to an 
unprecedented level of accuracy. Upon cooling, this new set of data unambiguously reveals a second 
order phase transition lowering the symmetry from tetragonal to orthorhombic at a temperature 
T c 2 = 136 K. The observation of an orthorhombic distortion of the a6-plane is attributed to nematic 
phase formation supported by local Jahn-Teller (JT) dynamical instability. At Tjv = 105 K, spins 
order and at T c i = 100 K the tetragonal structure is recovered with an elongated c-axis. 

PACS numbers: 78.20.-e, 78.20.Ls, 75.25.Dk, 71.70.Ej 


I. INTRODUCTION 

Si^VCb can be stabilized with vanadium ions in a 
4 + valence state in two different polytypes: (i) an or¬ 
thorhombic form whose peculiar magnetic dimer behav¬ 
ior has been previously investigated [1] and (ii) a tetrag¬ 
onal one that exhibits puzzling structural and magnetic 
transitions and is the object of the present study. The 
ground state of this Mott-Hubbard insulator has been 
studied experimentally [2-5] and theoretically [6-9] for 
a long time. Nevertheless, doubts remain about the na¬ 
ture of the order parameter. Upon cooling, a magneti¬ 
cally ordered ground state is achieved after a succession 
of structural and magnetic transitions, whose nature is 
not yet understood. Specific heat and magnetic suscepti¬ 
bility measurements have revealed the following magnetic 
and structural diagram: 

• Below 10 K the material exhibits a tetragonal struc¬ 
ture with a weak magnetic order [5, 6]. 

• Between 10 K and 100 K the system enters a differ¬ 
ent magnetic and orbital state, which is not ferro¬ 
magnetic and the nature of which has not yet been 
fully understood. 

• In the region between T c \ ~ 100 K and T c2 ~ 140 
K recent structural analysis suggests an orthorhom¬ 
bic distortion of the a6-plane [10] but the magnetic 
state remains not known. 

• Above T c 2 the material is again tetragonal with a 
slightly reduced c axis [3] and it exhibits a paramag¬ 
netic behaviour with a large Van Vleck contribution 
[9], 

The main candidates for the ground state that have 
been proposed until now are an octupolar ground state 
[8] and anti-ferro orbital ordering with muted magnetic 
moments due to spin-orbit coupling [4]. These two sce¬ 
narios are well supported by experiments for most of the 
temperature range between 10K and ambient tempera¬ 
ture. Only the splitting of the excited state measured 


with inelastic neutron scattering [11], while consistent 
with anti-ferro orbital ordering [4], is inconsistent with 
octupolar order[8]. However, neither of these two mod¬ 
els captures the two-stage ordering implied by the two 
peaks in the specific heat at 98 K and 127K. Two possible 
configurations account for the observed double structural 
transition at T c 1 and T c2 : (i) the coexistence of the low 
(LT) and the high temperature (HT) tetragonal phases, 
(ii) the onset of an intermediate orthorhombic distortion 
of the lattice associated to the transition from the HT to 
LT phase. 

This ambiguity was recently removed by Yamauchi et 
al. [10], who concluded, based on Le Bail fitting of X-ray 
powder diffraction data, that the intermediate phase is 
orthorhombic. They suggested the space group Immm 
as being the most likely, according to the group-subgroup 
relation to the IA/mmm space group of the tetragonal 
phase. Nevertheless, the structural model of the or¬ 
thorhombic phase was not resolved in Ref. 10, and higher 
resolution diffraction data are needed for that purpose. 

Here we present high resolution synchrotron X-ray 
diffraction data (XRD), and resolve the structural model 
of the orthorhombic phase that successfully fits the data, 
through a full pattern profile refinement performed using 
the Rietveld method. We also present theoretical argu¬ 
ments that the orbital ordering in Si^VC^is accompanied 
by an intermediate orthorhombic distortion stabilized by 
electron-lattice interactions. Examples of lattice medi¬ 
ated orbital ordering are well known in KCUF3 [12] or 
in LaMnC>3 [13, 14]. However, in those compounds e g 
electrons are involved, therefore the spin-orbit coupling 
does not play an essential role, contrary to the situation 

inSr 2 V0 4 . 


II. X-RAY POWDER DIFFRACTION DATA 

Synchrotron X-ray powder diffraction experiments 
were carried out at the high resolution diffraction beam¬ 
line Ill at the Diamond Light Source (details on experi¬ 
mental procedure and crystallographic data analysis are 



2 


given in appendix, details of the beamiline are reported 
in [15, 16]). The samples were prepared by reduction 
of Si' 4 V 20 g in a sealed quartz tube with zirconium as 
an oxygen getter [17]. No secondary phases nor traces of 
the orthorhombic phase were detected after 5 consecutive 
reduction cycles at 850 °C for 48h. 
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FIG. 1. a) Temperature evolution of selected Bragg reflec¬ 
tions. The two transition temperatures T c i and T c 2 are high¬ 
lighted with blue lines. Temperature intervals are 1 Kelvin 
from 20K to 40K, 5K from 45K to 85K, 1 Kelvin from 90K to 
139K, 10K from 140K to 190K and 20K from 200K to 280K. 

The temperature evolution of diffraction peaks in the 
aft-plane (200 in Fig. 1) and along the c-axis (002 and 004 
in Fig. 1) clearly indicates a long range aft-plane distor¬ 
tion through the splitting of the 200 reflection between 
T c 1 and T c 2 , and a first-order like c-axis expansion below 
T cl . During the Rietveld refinement (see details in ap¬ 
pendix) the symmetry of the structure was set to tetrag¬ 
onal in the high (above T c 2 ) and low (below T c i) temper¬ 
ature ranges which was sufficient to index all diffraction 
peaks. The diffraction pattern at 108 K, showing large 
departure from the two tetragonal phases (Fig. 2b) and 
slightly above TV = 105 K, was used to reliably identify 
the symmetry of the intermediate phase. A dichotomy in¬ 
dexing routine implemented in the program Fox [18] was 
used to index 20 low angle peaks belonging to Sr 2 V 04 . 


An orthorhombic cell corresponding to the deformation 
of the low temperature tetragonal cell (14/mmm) was 
immediately found, in agreement with a recent study in 
which long range orthorhombic distortion was also sug¬ 
gested [10]. The observed extinctions pointed only to a 
body centered lattice which is in agreement with 4 pos¬ 
sible space groups. The highest symmetry space group 
among them, /mram, was successfully used to solve the 
structure ab-initio by the global optimization method us¬ 
ing the program Fox , and modeling the structure by one 
Sr, one V and three O atoms. The resulting structural 
model corresponds to a deformation of the tetragonal 
structure with the VC>6 octahedron elongated along the 
c-axis, and slightly deformed in the equatorial plane (the 
parameters are provided in the appendix), which is in 
agreement with the loss of the 4-fold axis and resulting 
orthorhombic deformation. 

Upon cooling down to the transition at T c 2 , observed 
also in the entropy extracted from specific heat (see meth¬ 
ods in appendix and point 1 in Fig. 2d), the a and c 
lattice parameters contract (Fig. 2b). Upon lowering 
the temperature below T c 2 the a and ft-axis parameters 
separate, and the c-axis shrinks. The volume remains un¬ 
changed at the transition at T c2 (1% of volume collapse 
was reported in LaMnCU at the Jahn-Teller transition 
[13, 14, 19]), but it slightly expands (0.02%) upon cooling 
below T c i (Fig. 2a). The weak lattice contraction visible 
at the lowest temperatures is in agreement with the prox¬ 
imity of a transition to a magnetic ordered state [5]. The 
evolution of the full width at half maximum (FWHM) 
of selected diffraction peaks is shown in Fig. 2c. No 
significant broadening is observed at the orbital ordering 
temperature at T c 2 pointing toward a second order like 
transition. A discontinuity is visible around T c 2 in which 
the pattern was refined to two separate structural models 
(orthorhombic, red-green symbols, and tetragonal, blue 
symbols, in Fig. 2c, too close to each other to take into 
account a likely coexistence of both of them. 


III. DISCUSSION AND INTERPRETATION 

The XRD data indicates a local in-plane distortion 
which is not compatible with the orthorhombic distortion 
of the octahedra typically observed in KCuF 3 [12] and 
LaMnCU [19], where the 4 fold axis is preserved although 
each octahedron presents an orthorhombic distortion of 
its aft-plane. Indeed, the tiny aft splitting observed in 
our XRD measurements (Fig.4) is only compatible with 
a breaking of C4 symmetry, i.e. a long range orthorhom¬ 
bic distortion (see also ferro-distortion in Fig. 3d). 

We propose the following model: Due to spin-orbit 
coupling each V-ion has a spin-orbital state in which spin 
(s = 1/2) and orbital (effectively l = 1) moments com¬ 
bine to a j = 1/2 spin-orbital doublet with a muted mag¬ 
netic moment due to the factor of 2 difference in spin and 
orbital gyromagnetic ratio. Due to the absence of long- 
range order of the muted-moments above T c 2 , the crystal 
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FIG. 2. Temperature evolution of a) c/a,b ratio and volume 
of the cell, b) lattice parameters and c) Full Width Half Max¬ 
imum (FWHM) of selected diffraction peaks in the high tem¬ 
perature and low temperature tetragonal phases (blue sym¬ 
bols) and in the orthorhombic phase (red and green symbols). 
The two transition temperatures T c i and T c 2 are highlighted 
with dotted lines. Note that error bars in a) and b) are smaller 
than the symbols, d) entropy corrected for the phonon and 
thermal expension (see methods in appendix) 


structure is tetragonal. Below the orbital ordering tem¬ 
perature (T c2 ) coupling to Jahn-Teller (electron-lattice) 
modes (Fig. 3a) splits the adiabatic potential (Fig. 3b) 
driving the system to a local distortion. The locally dis¬ 
torted ions form a nematic phase. At the Neel tempera¬ 
ture T)v=105 K antiferromagnetic ordering sets in. The 
Jahn-Teller coupling switches off at T c i, thus allowing the 
system to come back to a bipartite square lattice. This 
restores the tetragonal symmetry of the lattice while the¬ 


oretically we expect a doubling of the unit cell. 







FIG. 3. a) representation of ionic displacement correspond¬ 
ing to the Q 2 , Qs and Q& Jahn-Teller distortion modes, b) 
Scheme of the adiabatic potential profile due to electron lat¬ 
tice interaction. Q is the distortion strength and Umin = 
ksToo the potential barrier. Scheme of c) antiferro- and d) 
ferro-distortions associated to coupling with each mode. 

In a previous publication [9], we have discussed the 
crystal field Hamiltonian including spin-orbit coupling 
and super-exchange interactions. The relevant orbital 
degrees of freedom is the set of degenerate vanadium 3d 
orbitals 3 d xz and 3 d y ~ for both spin-quantum numbers, 
which are mixed by spin-orbit interaction. These orbitals 
are coupled to the lattice degrees of freedom, Qj with 
j = 1,2,3,6 of which j = 2,3,6, depicted in Fig. 3a), 
are Jahn-Teller active. Only the Q 2 mode is sufficiently 
strongly coupled to be relevant in the present discussion. 
The adiabatic potential in reduced dimensionless units 
describing these couplings is 

u{q) = y - v(T, m) \/l + q 2 (1) 

where v(T,m) is the electron-lattice coupling constant 
normalized by the spin-orbit coupling constant, which 
varies as a function of temperature as well as the anti¬ 
ferromagnetic order parameter m (see Appendix). In 
particular v > 1 for T c 1 < T < T c 2 and v < 1 above and 
below this temperature range. On the level of the in¬ 
dividual V 4+ -ions, the displacements of the surrounding 
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FIG. 4. a) Experimental data of the Jahn-Teller order param¬ 
eter (a — b)/(a + 6 ), b) calculation of the coupled Jahn-Teller 
order parameter q and c) the anti-ferromagnetic order param¬ 
eter m using Eq. 7. The parameters are indicated in the figure 
and caracteristic temperatures T c \, Tn and T c 2 are marked as 
vertical dashed lines. 


oxygen ions are subject to zero-point motion. However, 
these motions are mechanically coupled between neigh¬ 
boring sites. The coordinate Q 2 describing the collective 
vibration of all ions then represents a macroscopic ob¬ 
servable, subject to the occurrence of spontaneous sym¬ 
metry breaking with two energetically equivalent minima 
for positive and negative Q 2 (Fig. 3b). The symmetry 
breaking disappears for high enough temperatures where 
thermal motion of the lattice overcomes the symmetry- 
broken state. The scenario is therefor as follows: For 
T > T c 2 the electrons occupy orbital momentum eigen¬ 
states d XZt(7 ± idyz tt7 , without long-range intersite mag¬ 
netic order. Below T c2 a collective Jahn-Teller effect 
breaks the electronic symmetry by shifting the orbital 
character towards d XZtC , or d yz>er while at the same time 
moving Q 2 away from the origin. A similar dynamical 
regime is present in LaMnCU but at much higher tem¬ 
perature ( Tjt = 700 K) due to strong electron-lattice 
coupling [13]. The static local Jahn-Teller effect usu¬ 
ally leads to an anti-ferro distorted order, which globally 
preserves the tetragonal (C4) symmetry (see Fig. 3c), 
but would induce a doubling of the unit cell in 06 -plane. 


This is not compatible with our experimental observa¬ 
tions where a long range distortion is observed (Fig. 2 
and Fig. 3d). We therefor conjecture that the dynamical 
Jahn-Teller local distortions enable a ferro-nematic long 
range order induced by interaction via the phonon field. 

The magnetic transition in S^VCU occurs near 100 K. 
Until now the experimental data have been interpreted 
as a first order transition, i.e. m(T ) jumps to a finite 
value at T cl (see Fig. 1 in ref.4). The spin-orbit coupling 
causes the orbitals to polarize as d xz ^ a ± id yzt<7 , and the 
Jahn-Teller instability is suppressed due to the fact that 
the parameter v jumps to a value smaller than 1 (see 
Appendix). Indeed the tetragonal symmetry is recovered 
below T c i as observed in the crystallographic data in Fig. 
2. However, the evolution of the a-axis and c-axis param¬ 
eters below T c i suggests another scenario where there is 
a second order anti-ferromagnetic transition at Tn- The 
magnetic order parameter m(T) than increases continu¬ 
ously (Fig. 4c). At T c i the condition v = 1 is reached, be¬ 
low which the Jahn-Teller distortion disappears. In this 
scenario there are three critical temperatures: T c i = 100 
K, Tn = 105 K and T c2 = Tjt = 136 K. The steps 
T c i K (point 1 in Fig. 2d) and T c2 (point 2 in Fig. 2d) 
in the entropy are associated with the tetragonal to or¬ 
thorhombic structural transitions that are visible in the 
crystallographic data around T c i K and T c2 . The temper¬ 
ature dependence of the Jahn-Teller distortion and the 
anti-ferromagnetic order are described by combining the 
Jahn-Teller adiabatic potential and a Ginzburg-Landau 
expression for the free energy of the anti-ferromagnetic 
order parameter, as detailed in the Appendix. The model 
calculation of the Jahn-Teller order parameter q (Fig. 4 
b) and the experimental a /6 axis ratio (Fig. 4 a) show 
qualitatively the same features: A gradual rise of q be¬ 
low T c 2 , sudden break at Tn and disappearance at T c i. 
These similarities suggest that the observed temperature 
dependence of the a /6 splitting is described by the com¬ 
bined effect of a Jahn-Teller distortion between T c i and 
T c 2 and the anti-ferromagnetic order of the spin-orbital 
mute moments below Tn- 

It is worth noting that a very similar scenario of 
interplay between nematic order and a small distortion 
of the crystal lattice was recently discussed in iron-based 
superconducting pnictides. The instability in pnictides is 
not related to JT dynamical effect but to Fermi surface 
nesting [20]. The resulting nematic phase induces a very 
similar long range orthorhombic distortion with a re¬ 
entrant tetragonal structure at low temperature [21]. A 
recent study on LaSrVCU reports that local Jahn-Teller 
distortions are responsible for a liquid orbital state [ 22 ]. 
This is indeed one of the signposts of the nematic order 
suggested by the present work. The splitting in the 
crystal field optical excitation spectrum at T c2 (Fig. 
5 in Ref. 4), originally attributed to the coexistence 
of two tetragonal phases, probably arises from the 
lifting of the degeneracy of the lowest Kramers doublet 
by the lowering of the symmetry below the JT transition. 
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IV. CONCLUSIONS 

In conclusion our high resolution powder diffraction 
measurements on tetragonal Sr 2 V 04 have revealed a 
second order transition to a weak orthorhombic distorted 
structure below T c 2 = 136 K. In the temperature window 
from T c 2 down to T c 1 = 100 K, superexchange and spin- 
orbit coupling, which tend to favor an anti-ferro-orbital 
ground state, compete with a nematic orbital ordering 
supported by the dynamical Jahn-Teller effect. Below 
the spin-orbital ordering temperature (Tn = 105 K), 
super-exchange and spin-orbit interaction screen the 
Jahn-Teller effect and at T c i the tetragonal symmetry is 
recovered. 
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VI. APPENDIX 
A. Experimental methods 
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FIG. 5. Rietveld refinement of the orthorhombic phase at 110 
K. Diffraction peaks for alluminium are marked with green 
stars. 

Structural studies were performed through X-ray 


parameter 

T=280 K 

a (A) 
c(A) 

3.83689(2) 

12.54294(12) 

R P 

0.164 

Rw P 

0.242 


parameter 

T=108 K 

a (A) 

3.82642(3) 

b(A) 

3.83438(3) 

c (A) 

12.52961(8) 

R P 

0.155 

Rw P 

0.199 


parameter 

T=20 I< 

a (A) 

c (A) 

3.82688(3) 

12.54771(16) 

R P 

0.277 

Rwp 

0.368 


TABLE I. lattice parameters and R-factors for the 3 crystal¬ 
lographic structures at 280K (high temperature (HT) tetrag¬ 
onal phase), 108 K (distorded orthorhombic phase), 20K (low 
temperature (LT) tetragonal phase). 


diffraction experiments on powder samples. Measure¬ 
ments were done at synchrotron facility at Diamond 
Light Source, by using the Ill beamline [15]. The X- 
ray beam was monochromatized at E = 15 keV (A = 
0.82713 A). The high-resolution detector was made of five 
multiple-analyzing-crystals and spanned a 26 range from 
0.02° to 150.00° with a Ad step of 0.004° [16]. The sam¬ 
ple was made of a finely ground Sr 2 V 04 powder, homoge¬ 
neously coating a thin A1 wire perpendicular to the beam 
direction and the scattering vector. Thanks to this sam¬ 
ple holder, any preferred orientation of Sr 2 V 04 grains was 
checked to be absent. The sample holder was mounted in¬ 
side a cryogen-free refrigerator and the temperature was 
controlled in the range from 20 to 300 K, while XRD 
diffraction was acquired at constant temperature. Small 
temperature steps (IK) were chosen over the range where 
structural transitions were expected to occur (20 — 40 K 
and 90 — 136 K). With the chosen number of T steps, the 
acquisition time at each T and the heating rate (0.1 — 1 
K/min), the whole diffraction experiment took about 25 
hours. During the allocated beamtime, we could repeat 
the experiment twice, on two different powder samples, 
for checking the reproducibility. 

The diffraction data analysis must carefully take into 
account the contribution of the sample holder which gives 
the strongest contribution to scattering. Nevertheless, 
Fig. 5 shows a very strong signal to noise ratio of low- 
intensity diffraction peaks at large angle, which proves 
the high quality of the experimental data and allows ex¬ 
tracting reliable information. 

A full pattern profile refinement based on the Rietveld 
method was used by means of the FullProf Suite program 
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2.05 [23]. The first parameters to be refined were those 
concerning the A1 sample holder. As a matter of fact, 
the A1 wire was found to have a filament-like (extrusion) 
texture, and therefore to exhibit a strong preferred ori¬ 
entation. 

Three known phases contribute to the diffraction pat¬ 
tern and were included in the pattern profile refinement: 
A1 (the sample holder), Sr 2 V 04 (either tetragonal or or¬ 
thorhombic depending on the temperature range, as ex¬ 
plained in the following), and Sr 3 V 2 0 8 (the only possible 
impurity phase present in the starting powder). 26 pa¬ 
rameters were refined in total, namely the scale factors, 
lattice parameters, peak shape (pseudo-Voigt) for A1 and 
Sr 2 VC> 4 , preferred orientation and asymmetry correction 
for A1 only, and zeroshift. After having refined the lat¬ 
tice parameters of all phases, the coordinate of the apical 
oxygen of the VC>6 octahedron was refined as well. The 
background was drawn by linear interpolation between 
manually selected points. 

The profile refinement strategy was the following: in a 
first step we refined the scale, lattice, texture and asym¬ 
metry parameters of Al, and the zeroshift. Then, with 
all these parameters constrained to their convergence val¬ 
ues, we refined scale and lattice parameters of the oxides. 
In a third step we refined lattice and shape parameters 
of Sr 2 VC> 40 nly, with all the others constrained at their 
primary convergence values, before a last step in which 
all parameters were free and refined together to the final 
convergence values. 

The R factors are reported to be in the range 0.15 - 
0.4 (see Table I) and could seem to be rather high at 
first glance. The quality of the fits is mainly limited 
by i) the strong filament-like texture of the Al sample 
holder; ii) the impossibility to refine two (likely) coex¬ 
isting tetragonal and orthorhombic Sr 2 V 04 phases near 
to the transitions temperatures; iii) the fixed composi¬ 
tion, and therefore the site occupation and position, of all 
phases; iv) the fixed zeroshift parameter over the whole 
temperature range. Despite this affects the R-factor val¬ 
ues, it has no influence on the lattice constants extracted 
from the structural refinements. 

The lattice parameters of Al were found to be in excel¬ 
lent agreement with the expected calibration curve. The 
secondary phase, Sr 3 V 2 0s, was found to be present at 
1.5 — 2% volume fraction, in agreement with the pre¬ 
liminary laboratory XRD experiments that could not 
identify any secondary phases within the limit of in¬ 
strument accuracy. Three different temperature ranges 
were chosen, in which three different structure models 
for Sr 2 V 04 were used for the pattern refinement. Below 
T c i, and above T c2 , a tetragonal structure with space 
group 74 /rmrirri was chosen, whereas an orthorhombic 
one with a space group Immm was used in the range 
T c i~T C 2 . The orthorhombic distortion at intermediate 
temperature was obtained from ab-initio structure so¬ 
lution. At crossover temperatures, corresponding to the 
structural transitions, each Sr 2 V 04 phase was kept well 
beyond its stability limit in order to check the continu¬ 


ous overlap of the lattice parameters. Correspondingly, 
a divergence of the peak shape parameters (determining 
the full width at half maximum, FWHM) of the tetrago¬ 
nal and orthorhombic phases was found, as expected (see 
Fig. 2 c). 

A summary of the lattice parameters is given in Ta¬ 
ble I for the 3 crystallographic structures at 280K (High 
temperature (HT) tetragonal phase), 105 K (most dis¬ 
torted orthorhombic phase), 20K (Low temperature (LT) 
tetragonal phase) 

The entropy shown in Fig. 2 d is the integration of 
residual specific heat after subtraction of the phonon con¬ 
tributions which we modeled with one Einstein and one 
Debye mode, and a linear term representing thermal ex¬ 
pansion. 


B. Details of the model 


Due to the tetragonal crystal field the Hilbert 
space relevant for the low energy properties is 
spanned by the degenerate set of single-electron states 
{\yz, +), \yz, —), \xz, +), \xz, —)}. Due to the symmetry 
properties of this manifold it only couples to Q 2 and Qq. 
We have verified numerically that the coupling to Qq is 
too small to be relevant to the present discussion. In 
the sequel we therefor ignore the coupling to Qq. The 
effect of spin-orbit coupling is to mix these states and 
form a set of non-degenerate spin-orbital doublets. The 
electron-lattice interaction and the spin-orbit coupling at 
a given site j are described by the Hamiltonian [9, 24] 


2 M + 2 


V 

~2 i 


P + + E 


Q + A l z s z 


( 2 ) 


Here k is the vibrational force constant, V is the electron- 
lattice coupling constants and P the momentum of the 
vibrational coordinate. Expressing the electronic hamil- 
tonian on the basis d± i jCr = d XZt<7 ± id yZt(7 we obtain a 
blockwise diagonal matrix with 2x2 blocks 


p2 

2 M 


kQ 2 
2 


( A/2 iVQ \ 
{-iVQ -A/2 ) ' 


( 3 ) 


We have to take into account the inter-site coupling of the 
Jahn-Teller active sites. For the purpose of the present 
discussion we assume a strong coupling, implying uni¬ 
form displacement of all Q(j) corresponding to the dif¬ 
ferent lattice sites j. We define the dimensionless pa¬ 
rameters v = 2V 2 /(kX), h = 477V 2 /(fcA 2 ), q = 2QV/X, 
p = PX/(2HV), p = MkX 4 /(16H 2 V 4 ), and obtain 


h = 


2/i 


r 

2 



( 4 ) 


Straightforward diagonalization, and treating q as a clas¬ 
sical variable, constitutes the adiabatic potential 


u{q) 



V 1 + q 2 


( 5 ) 
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The first important observation is, that for v < 1 the 
absolute minimum occurs at q = 0, hence symmetry 
breaking requires that v > 1. In the present mate¬ 
rial the spin-orbit parameter is A = —30 meV[9]. Fur¬ 
thermore, using experimentally-based data for RbMnF 3 
[25]) and KMgF 3 [26, 27], we'obtain V = 3 1 / 2 3 V E /2 = 
—0.26 ± 0.05 eV/A, and the relevant force constant is 
k = 2eV/A 2 . With these parameters the coupling con¬ 
stant v = 1.5 ± 0.5, which is above the critical value 
v = 1. It is furthermore natural to suppose that v has 
a temperature dependence due to thermal fluctuations, 
which limits the range where v(T) > 1. 

Anti-ferro-magnetism orders the spins along the the c- 
axis. Due to the spin-orbit coupling this effectively sup¬ 
presses the Jahn-Teller coupling. The minimal model 
describing this state of affairs is 

v{T,m) = 1+0(1 -am 2 )(l-T/T c2 ) (6) 

where a is a constant describing the coupling between the 
Jahn-Teller order parameters, and the anti-ferromagnetic 
order parameter m. The coupled Jahn-Teller and anti- 
ferromagnetic order are described by the free energy 

q 2 

u(q, m) = — - v(T,m)^/l + q 2 + f{m) 
f(m) = f 2 m 2 + / 4 m 4 5 + / 6 m 6 7 8 9 10 (7) 


The free energy minimum corresponds to the condi¬ 
tion du/dq = du/dm = 0. The Jahn-Teller transition 
at T c2 (where in the case of Sr 2 V 04 there is no anti¬ 
ferromagnetic order) follows from the condition 

m(T c2 ) = 0 

v(T c2 , 0) = 1 ( 8 ) 

A second order magnetic transition is described by the set 
of parameters fe = 0 , / 4 = 1/2 and f 2 = q(T /T m — 1). 
The anti-ferromagnetic transition occurs at Tjv, where 
T)v < T m due to the competition between Jahn-Teller 
ordering and anti-ferromagnetism. Lowering the temper¬ 
ature below Tjv, the anti-ferromagnetism progressively 
suppresses the Jan-Teller ordering. This process is com¬ 
pleted at the temperature T cl , which follows from the 
self-consistent relation 

v(T cl ,m(T cl )) = l (9) 

For illustration we show in Fig. 4 the coupled order pa¬ 
rameters q and m as a function of temperature using Eq. 
7. The input parameters a = 0.35, /3 = 1, 7 = 50, 
T c 2 = 136 K, T m = 105 K, were chosen such as to ob¬ 
tain values for the transition temperatures of the coupled 
system Tjv (105 K), and T c 1 (100 K) close to the experi¬ 
mental ones. 


[1] J. Deisenhofer, S. Schaile, J. Teyssier, Z. Wang, M. Hem- 
mida, H.-A. K. von Nidda, R. M. Eremina, M. V. Eremin, 
R. Viennois, E. Giannini, D. van der Marel, and A. Loidl, 
Physical Review B 86, 214417 (2012). 

[2] M. Cyrot, B. Lambert-Andron, J. Soubeyroux, M. Roy, 
P. Dehauht, F. Cyrot-Lackmann, G. Fourcaudot, 
J. Beille, and J. Tholence, Journal of Solid State Chem¬ 
istry 85, 321 (1990). 

[3] H. D. Zhou, B. S. Conner, L. Balicas, and C. R. Wiebe, 
Physical Review Letters 99, 136403 (2007). 

[4] J. Teyssier, R. Viennois, E. Giannini, R. M. Erem¬ 
ina, A. Gunther, J. Deisenhofer, M. V. Eremin, and 
D. van der Marel, Physical Review B 84, 205130 (2011). 

[5] J. Sugiyama, H. Nozaki, I. Umegaki, W. Higemoto, E. J. 
Ansaldo, J. H. Brewer, H. Sakurai, T.-H. Kao, H.-D. 
Yang, and M. Mansson, Physical Review B 89, 020402 
(2014). 

[6] W. Pickett, D. Singh, D. Papaconstantopoulos, 
H. Krakauer, M. Cyrot, and F. Cyrot-Lackmann, Phys- 
ica C: Superconductivity 162—164, Part 2, 1433 (1989). 

[7] Y. Imai, I. Solovyev, and M. Imada, Physical Review 
Letters 95, 176405 (2005). 

[8] G. Jackeli and G. Khaliullin, Physical Review Letters 
103, 067205 (2009). 

[9] M. V. Eremin, J. Deisenhofer, R. M. Eremina, 
J. Teyssier, D. van der Marel, and A. Loidl, Physical 
Review B 84, 212407 (2011). 

[10] I. Yamauchi, K. Nawa, M. Hiraishi, M. Miyazaki, 
A. Koda, K. M. Kojima, R. Kadono, H. Nakao, R. Ku- 
mai, Y. Murakami, H. Ueda, K. Yoshimura, and M. Taki- 
gawa, Phys. Rev. B 92, 064408 (2015). 


[11] H. D. Zhou, Y. J. Jo, J. Fiore Carpino, G. J. Munoz, 
C. R. Wiebe, J. G. Cheng, F. Rivadulla, and D. T. 
Adroja, Phys. Rev. B 81, 212401 (2010). 

[12] J. C. T. Lee, S. Yuan, S. Lai, Y. I. Joe, Y. Gan, 
S. Smadici, K. Finkelstein, Y. Feng, A. Rusydi, P. M. 
Goldbart, S. Lance Cooper, and P. Abbamonte, Nature 
Physics 8, 63 (2012). 

[13] M. R. Ahmed and G. A. Gehring, Physical Review B 79, 
174106 (2009). 

[14] A. Trokiner, S. Verkhovskii, A. Gerashenko, Z. Volkova, 
O. Anikeenok, K. Mikhalev, M. Eremin, and L. Pinsard- 
Gaudart, Physical Review B 87, 125142 (2013). 

[15] J. E. Parker, S. P. Thompson, T. M. Cobb, F. Yuan, 
J. Potter, A. R. Lennie, S. Alexander, C. J. Tighe, J. A. 
Darr, J. C. Cockcroft, and C. C. Tang, Journal of Ap¬ 
plied Crystallography 44, 102 (2011). 

[16] S. P. Thompson, J. E. Parker, J. Potter, T. P. Hill, 

A. Birt, T. M. Cobb, F. Yuan, and C. C. Tang, Review 
of scientific instruments 80 (2009), 10.1063/1.3167217. 

[17] R. Viennois, E. Giannini, J. Teyssier, J. Elia, J. Deisen¬ 
hofer, and D. V. der Marel, Journal of Physics: Confer¬ 
ence Series 200, 012219 (2010). 

[18] V. Favre-Nicolin and R. Cerny, Journal of Applied Crys¬ 
tallography 35, 734 (2002). 

[19] T. Chatterji, F. Fauth, B. Ouladdiaf, P. Mandal, and 

B. Ghosh, Physical Review B 68, 052406 (2003). 

[20] R. M. Fernandes, A. V. Chubukov, and J. Schmalian, 
Nature Physics 10, 97 (2014). 

[21] S. Avci, O. Chmaissem, J. M. Allred, S. Rosenkranz, 
I. Eremin, A. V. Chubukov, D. E. Bugaris, D. Y. Chung, 
M. G. Kanatzidis, J.-P. Castellan, J. A. Schlueter, 



H. Claus, D. D. Khalyavin, P. Manuel, A. Daoud- 
Aladine, and R. Osborn, Nature Communications 5 
(2014). 

[22] Z. L. Dun, V. O. Garlea, C. Yu, Y. Ren, E. S. Choi, H. M. 
Zhang, S. Dong, and H. D. Zhou, Physical Review B 89, 
235131 (2014). 

[23] J. Rodrguez-Carvajal, Physica B: Condensed Matter 
192, 55 (1993). 


[24] A. Abragam and B. Bleaney, Electron paramagnetic res¬ 
onance of transition ions (Clarendon P., 1970). 

[25] E. I. Solomon and D. S. McClure, Physical Review B 9, 
4690 (1974). 

[26] S. Muramatsu and N. Sakamoto, Journal of the Physical 
Society of Japan 46, 1273 (1979). 

[27] M. G. Brik, N. M. Avram, and I. Tanaka, physica status 
solidi (b) 241, 2982 (2004). 



